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Recent experiments performed at high pressures indicate 
that ruthenium can support unusually high concentrations of 
oxygen at the surface. To investigate the structure and stabil- 
ity of high coverage oxygen structures, we performed density 
functional theory calculations, within the generalized gradient 
approximation, for O adlayers on Ru (0001) from low cover- 
age up to a full monolayer. We achieve quantitative agree- 
ment with previous low energy electron diffraction intensity 
analyses for the (2 x 2) and (2 x 1) phases and predict that 
an O adlayer with a (1 x 1) periodicity and coverage 0=1 can 
form on Ru(OOOl), where the O adatoms occupy hep- hollow 
sites. 

PACS numbers: 68.35.-p, 82.65.My 



I. INTRODUCTION 

The interaction of oxygen with metal surfaces forms 
the basis of many technologically important processes, for 
example, bulk oxidation, corrosion, and heterogeneous 
catalysis. It is therefore of great interest to obtain a 
detailed understanding of the changes in the atomic and 
electronic structure that oxygen adsorption often induces 
due to the formation of strong chemical bondsElfl. The be- 
havior of O on metal surfaces is quite varied and depends 
markedly on the coverage and temperature, and on the 
orientation of the surface of the particular metal. Gen- 
erally, the close-packed surfaces are more stable against 
reconstruction; often, however, significant atomic eelax- 
ations of the substrate are induced by O adsorptionEl. On 
Zr, the O atoms apparently form an underlayerEI even for 
the close-packed (0001) surface. At higher coverages of 
oxygen, at elevated temperatures, oxick-like structures 
can form on a number of metal surfacesB. 

From recent experiments of the catalytic oxidation of 
carbon monoxide, performed at high pressures, there is 
evidence that Ru (0001) can suppoil Jjnusually high con- 
centrations of oxygen at the surfacaj'Q. In order to inves- 
tigate the structure and stability of high coverage oxy- 
gen structures, we performed extensive density functional 
theory calculations for various O adlayers on Ru(W)Ol). 
In parijirular, for the two ordered phases, (2 x 2)0 and 
(2 X 1)13, which form at room temperature under ultra 
high vacuum (UHV) conditions for coverages = 1/4 
and O = 1/2, respectively, as well as for an artificial 
(2 X 2)-30 adlayer with coverage = 3/4, and for sev- 
eral higher coverage (1 x l)-0 structures with = 1. 
Here, O is defined to be the ratio of the number of adsor- 
bate atoms to the number of atoms in an ideal substrate 



layer. The calculations reveal that although a (1 x 1) 
phase is not observed to form under UHV conditions us- 
ing molecular oxygen, perhaps due to the presence of 
activation energy barriers for dissociation of O2, the ad- 
sorption ofOina(lxl) adlayer structure with coverage 
O = 1 is exothermic and should be able to form when 
energy barriers can be overcome. 

II. CALCULATION METHOD 

We use density functional theory (DFT) and employ 
the generalized gradient approximation (GGA) of Perdew 
et alB for the exchange-correlation functional. The sur- 
face is modelled using the supercell approach and we em- 
ploy ab initio, fully separable pseudouptentials created by 
the scheme of Troullier and Martingj, and a plane- wave 
basis set for expansion of the electronic wave functions. 
The calculations are performed using a (2 x 2) surface 
unit cell and four layers of Ru (0001) with a vacuum re- 
gion corresponding to thirteen such-layers. The O atoms 
are adsorbed on one side of the slablij. The energy cut-off 
is taken to be 40 Ry aiith three special k-points uj-the 
surface Brillouin zondiil. The calculation schemeE£l af- 
fords simultaneous relaxation of the electrons and atoms 
using damped Newton dynamics. We relax the positions 
of the O atoms and the Ru atoms in the top two layers, 
keeping the lower two Ru layers fixed. For the higher 
coverage (1 x l)-0 structures we use a larger energy cut- 
off of 60 Ry and 14 k-points in the irreducible part of the 
surface Brillouin zone of a (1 x 1) surface unit cell. 

III. O ON RU (0001) 

A. (2 x 2)-0/Ru (0001) 

We performed calculations for O in the fcc-hoUow site 
(no atom in the layer beneath the site) and in the hcp- 
hoUow site. From our calculations we find that the hcp- 
hoUow site is energetically clearly favorable. This is 
in agreement with the site determined by a dynamical 
LEED intensity analysiaj. The binding energy of O (rel- 
ative to a free O atom) in the hcp-hoUow site is 5.55 eV 
and in the fcc-hollow site it is 5.12 eV (see Tab. I). We 
can compare this result with calculations for a single O 
adatom op_an eighteen atom cluster model representing 
Ru (0001)113, using an atom superposition and electron 
delocalization molecular orbital ASED-MO method. In 
this study it was found that the fcc-hollow site is the most 
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favorable with a binding energy of 5.6 eV. The hcp-hoUow 
site was the next most favorable with a binding energy 
of 5.3 eV. Such cluster geometries are often regarded as 
representing the situation of an isolated adsorbate on an 
extended surface. Our calculated coverage dependence of 
the adsorption site of O, discussed below (see Fig. 4a), 
however, suggests that an isolated O adatom will occupy 
the hcp-hoUow site. 

The atomic geometry of (2 x 2)-0/Ru (0001) is dis- 
played in Fig. 1. In Tab. II we compare the theoreti- 
cally obtained structutal parameters with those obtained 
by the LEED analysisO. For ease of comparison, we use 
the same notation as in the LEED study. We note that 
the theoretical lattice constant is 1.9 % larger than the 
experimental one. Rather than giving a detailed com- 
parison of all the substrate relaxations, we simply refer 
to Tab. II, from which the high level of agreement with 
the LEED analysis can immediately be assessed. The 
theoretically obtained 0-Ru bond length of 2.08 A is 
slightly longer than the LEED-determined value of 2.03 
A . The first Ru-Ru interlayer spacing is found to be con- 
tracted by 2.7 % with respect to the bulk value (using 
the centers of gravity of the first and second buckled Ru 
layers) . This agrees well with the value determined from 
the LEED analysis of 2.1 %. Such a contraction is in 
contrast to many similar systems formed by adsorption 
of strongly electronegative species, where instead there is 
often an expansion induced by the adsorbatetl. We find 
the top Ru-Ru interlayer spacing of Ru (0001) to be con- 
tracted by 2.5 p% which is close to the LEED-determined 
value of 2.3 The contraction of the top interlayer 

spacing of the clean surface is therefore not removed by 
oxygen adsorption at 1/4 monolayer. 

B. (2 X l)-O/Ru(0001) 

At higher oxygen coverage, namely, 8 = 1/2, a(2x2) 
LEED pattern is observed which corresponds to three 
rotated domains, each of (2 x 1) periodicityQ. We per- 
formed calculations for O in the fee- and hcp-hoUow sites. 
The hcp-hoUow site is again found to be energetically the 
most favorable with a binding energy of 5.28 eV; that for 
the fcc-hoUow site is 5.00 eV (see Tab. I). The theoretical 
identification of the hep- hollow site for O in the (2x1) 
structure is in accerd with the LEED determination for 
the adsorption siteQ. 

The atomic structure of (2 x l)-0/Ru (0001) is depicted 
in Fig. 2. The O atoms adsorb in "off" hcp-hoUow sites, 
i.e., they are displaced from the center of the hcp-hollow 
site. In addition, complex O induced relaxations of the 
substrate occur, including "row-pairing" and buckling of 
the substrate layers. The determined 0-Ru bond length, 
and the lateral and vertical relaxations are given in Tab. 
Ill where they are compared with the results obtained 
from the LEED analysisQ. Again, it can quickly be seen 
that quantitative agreement is achieved. We do note one 



deviation however: The directions of the lateral displace- 
ment of atom D, A(i||(D), have the opposite signs. That 
is, we obtain row-pairing of the Ru atoms in the second 
Ru layer, as well as in the first layer, and the LEED anal- 
ysis does not. We found that relaxing the third Ru layer 
does not change this result. The 0-Ru bond length of 
2.07 A is very similar to that which we determined for 
the lower coverage (2 x 2) structure which was 2.08 A. 
The value is again somewhat larger than that of 2.02 A as 
obtained from the LEED analysis. The first two Ru-Ru 
interlayer spacings, defined with respect to the centers of 
gravity of the buckled atomic layers, correspond to the 
bulk value to within 0.01 A , for both the DFT-GGA and 
LEED results. 



C. (1 X l)-0/Ru (0001) 

We now investigate the structure and stability of 
(1 X l)-0 adlayers with coverage = 1. We performed 
calculations for O in the on-top, bridge, fee-, and hcp- 
hollow sites. The obtained binding energies are collected 
in Tab. I. The hcp-hollow site is energetically preferred 
with a binding energy of 4.87 eV. The fcc-hoUow site 
is the next most favorable having a binding energy of 
4.81 eV. Table I also lists the binding energy differences 
for the adsorption sites tested, with respect to the bind- 
ing energy of O in the respective hcp-hollow sites. The 
atomic geometry of (1 x l)-0 /Ru (0001) is shown in Fig. 3 
and the 0-Ru bond length and structural parameters are 
given in Tab. IV. It is noticeable that the 0-Ru bond 
length of 2.03 A is slightly shorter than that of the lower 
coverage structures. The first Ru-Ru interlayer spacing 
is found to be expanded by 2.7 %. 

For comparison we also performed calculations for O 
in the hcp-hollow site using an energy cut-off ofj-AO Ry 
and three special k-points in the Brillouin zoneEJ. As 
can be seen from Tab. IV, the resulting binding energy 
is 4.84 eV, i.e., only 0.03 eV less than the result of the 
more accurate calculation, and the resulting structural 
parameters differ by at most 0.02 A . 

The value of the binding energy of O in the hcp-hollow 
site on Ru (0001) at coverage 9 = 1 shows that the ad- 
sorption is exothermic and indicates that the (1x1) 
adlayer structure should be able to form. That is, the 
binding energy is larger (by « 1.8 eV per atom) than 
that which the O atoms have in O2. Under UHV con- 
ditions, however, the (2 x 1) phase is the terminal one. 
The reason that the (1 x 1) structure does not form under 
UHV conditions could be due to a kinetic hindering of the 
dissociation of O2 induced by the (2 x 1) structure. Inter- 
estingly, on a stepped Ru (0001) surface, the formation of 
a (1 X 1) structure for covecage = 1 has been reported, 
which is stable to 600 KEj. On the stepped surface, it 
is possible that step edges may act as sites over which 
dissociation of O2 can occur. If atomic, as opposed to 
molecular, oxygen would be used then it is probable that 
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the (1 X 1) phase will also be observed on Ru (0001) under 
UHV. This result could have implications for heteroge- 
neous catalytic reactions in which dissociative adsorption 
of O2 is a necessary reaction step (often rate-limiting) in 
that if atomic oxygen would be used the kinetics may be 
greatly altered; it also raises the question if other metal 
surfaces can also support high O coverages. 

IV. COVERAGE DEPENDENCE 

As we have seen from above, the hcp-hollow site is the 
preferred adsorption site for O on Ru (0001) at all the 
coverages investigated. This is consistent with the com- 
mon trend that atoms strongly chemisorbed on transition 
metal surfaces usually occupy the site that the next sub- 
strate layer would occupy. 

In Fig. 4a the binding energy of O in the fee- and hcp- 
hollow sites as a function of coverage is displayed. With 
increasing coverage the binding energy becomes less fa- 
vorable which reflects a repulsive interaction between the 
adsorbates and implies that no island formation is ex- 
pected to occur in the the coverage regime of 9 = 1/4 
to Q = 1. Figure 4a also shows that the difference in 
binding energy between the fee- and hcp-hollow sites be- 
comes less, as the coverage increases. As noted above, we 
also performed calculations for a structure with coverage 
= 3/4. In this structure we placed O atoms in hcp- 
hollow sites in the (2 x 2) surface unit cell. The 0-Ru 
bond length is 2.07 A which is similar to that of the two 
lower coverage structures, and the first Ru-Ru interlayer 
spacing is expanded by 1.8 %. The corresponding work 
function change as a function of coverage is shown in 
Fig. 4b. The experimental results of Surnev et allB are 
included for comparison where good agreement between 
theory and experiment is obtained. The large increase in 
work function reflects the strong electronegative nature 
of O. 



V. CONCLUSION 

From our calculations we achieve quantitative agree- 
ment with previous structural determinations for (2 x 2)- 
O and (2 x l)-0/Ru (0001) structures, with respect to 
the preferred adsorption site and the detailed atomic po- 
sitions in the surface region. One minor difference is 
found for the (2x1) structure which would be worth a 
re-analysis of the LEED data. We predict that an O ad- 
layer with a (1 X 1) periodicity and coverage Q=l should 
be able to form on Ru(OOOl), where the O atoms oc- 
cupy hcp-hollow sites. For the coverage regime 8=1/4 
to 0=1 there is no indication of island formation. The 
corresponding work function change agrees well with ex- 
periment. 
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FIG. 1. Top view (a) and side view (b) of the atomic geom- 
etry of (2 X 2)-0/Ru (0001). The arrows indicate the direction 
of the displacements of the substrate atoms with respect to 
the bulk terminated positions. The dashed line in (a) indi- 
cates the plane of the cross-section used in (b). Small dark 
grey circles represent oxygen atoms and large grey circles rep- 
resent Ru atoms. Interlayer spacings are given in Angstrom. 

FIG. 2. Top view (a) and side view (b) of the atomic geom- 
etry of (2 X l)-0/Ru (0001). The arrows indicate the direction 
of the atomic displacements. The dashed line in (a) indicates 
the plane of the cross-section used in (b). Small dark grey cir- 
cles represent oxygen atoms and large grey circles represent 
Ru atoms. Interlayer spacings are given in Angstrom. 
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FIG. 3. Top view (a) and side view (b) of the atomic geom- 
etry of (1 X l)-0/Ru (0001) witli O in the hcp-hoUow site. The 
arrows indicate the direction of the displacements of the sub- 
strate atoms with respect to the bulk terminated positions. 
Small dark grey circles represent oxygen atoms and large grey 
circles represent Ru atoms. Interlayer spacings are given in 
Angstrom. 

FIG. 4. (a) Binding energy of O on Ru (0001) for O in the 
fee- (dashed hue) and in the hep- hollow site (continuous line) . 
(b) Work function change as a function of coverage, Q. The 
experimental results, shown as open circles, are from Ref.LJ. 



bites tested 


on-top 


bridge 


icc-hollow 


hcp-nollow 


(2 X 2) 






5.12 


5.55 


AE 






0.43 


0.00 


(2x1) 






5.00 


5.28 


AE 






0.28 


0.00 


(ixl) 


3.62 


3.93 


4.81 


4.87 


AE 


1.25 


0.94 


0.06 


0.00 



TABLE I. Binding energies (in eV) of O on Ru (0001) relative to the free O atom for the surface structures investigated. 
The binding energy differences, A E, defined relative to the value for the respective hcp-hollow sites, are also given. 





(2 X 2)-0/Ru (0001) 


Structural parameters (A) 


O-Ru 


Ad||(A) 


Ad||(D) 


Ad. (A) 


Ad. (B) 


Ad. (C) 


LEED 


2.03 


0.09 


0.01 


-0.05 


-0.12 


-0.08 


DFT-GGA 


2.08 


0.07 


0.01 


-0.04 


-0.11 


-0.03 



TABLE 11. Structural parameters for (2 x 2)-0/Ru (0001) with O in the hcp-hollow site. A(i|| and Adz represent, respectively, 
lateral and vertical relaxations with respect to the bulk positions of the atoms indicated by the letters in parenthesis (see Fig. 1). 
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('9 V 1 
(Z X 1) 


-O/Ru (0001) 




Structural parameters 


(A) 


0-Ru 


Adz(A) 


Ad.(B) 


Ad.(C) 


Ad.(D) 


LEED 




2.02 


-0.03 


0.04 


-0.01 


0.02 


DFT-GGA 




2.07 


-0.06 


0.01 


-0.02 


0.00 


Structural parameters 


(A) 


Ad||(0) 


Ad||(A) 


Ad||(B) 


Ad|i(C) 


Ad||(D) 


LEED 




-0.06 


-0.07 


0.05 


0.05 


0.04 


DFT-GGA 




-0.03 


-0.02 


0.05 


0.01 


-0.02 



TABLE III. Structural parameters for (2x l)-0/Ru (0001) with O in the hcp-hoUow site. Ad|| and Adz represent, respectively, 
lateral and vertical relaxations with respect to the bulk positions of the atoms indicated by the letters in parenthesis (see Fig. 2). 









(i 


X Ij-O/Hii 


(DODi) 




Structural parameters (A) 


O-Ru 


dz,i 


dz,2 


dz,3 


dz.bulk 


Ei, (eV) 


DFT-GGA (40 Ry) 


2.04 


1.28 


2.27 


2.19 


2.19 


4.84 


DFT-GGA (60 Ry) 


2.03 


1.26 


2.25 


2.17 


2.19 


4.87 



TABLE IV. Structural parameters for (1 x l)-0/Ru (0001) with O in the hcp-hoUow site. 0-Ru, dz, and E;, represent, the 
O-Ru bond length, the interlayer spacings, and binding energy, respectively. 
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12 




0.09 
0.01 
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> 

(D 

0) 

bX) 

• i-H 

• i-H 

PQ 



4.6 p 

4.8 

5.0 

5.2 

5.4 

5.6 

5.8 



(a) 



0.0 0.2 0.4 0.6 0.8 
Coverage 



1.0 




0.2 0.4 0.6 0.8 
Coverage 
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